Equatoiral winds on Saturn and the Stratospheric Oscillation 
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The zonal jets on the giant planets are generally thought to be stable with time 1 \ Recently, 
there are still some debates about the general thought 4 ' 5 . Here, we report a significant 
temporal variation of the equatorial jet at high-altitude on Saturn. Long-term (2004-2009) 
observations by Cassini reveal that wind speed at the 60-mbar level increased from 270 m/s 
in 2004 to 290 m/s in 2008, while the wind speed has been mostly constant over time at the 
500-mbar level in the southern equatorial region. The Cassini observations further reveal 
that the equatorial jet intensified ~ 60 m/s in the stratosphere (1-5 mbar) from 2005 to 
2008. The fact that the wind acceleration is weaker at the 60-mbar level (~ 20 m/s) than at 
the I -mbar level (~ 60 m/s) demonstrates that the equatorial oscillation” 7 is damped when it 
propagates downwards to the tropopause around 60 mbar. The direct measurement of the 
varying equatorial jet around the tropopause also serves as a key boundary condition when 
deriving the thermal wind fields in the stratosphere. 

Does the equatorial jet vary with time on the giant planets? A comparison 4 of zonal winds 
between Voyager (1981-1982) and Hubble Space Telescope (1995-2002) observations suggested 
that velocity of Saturn’s equatorial jet has decreased with time. Conversely, an analysis 5 of 
observations by the Imaging Science Subsystem (ISS) on Cassini argued that the decrease of the 
equatorial jet discussed in the previous study 3 could instead be explained by changes in altitude 
of the tracked clouds within a zonal wind structure with vertical shear. Studies 8,9 based on data 
from the Composite Infrared Spectrometer (CIRS) on Cassini determine that vertical shear is 
present in the equatorial zonal winds. However, other studies 10,11 suggest that the vertical shear is 
not strong enough to explain the wind decrease. The above debates reveal the need to further 
characterize the equatorial jet on Saturn with altitude and time. Such results also are critical for 
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discriminating between different theories of large-scale circulation on the giant planets because 
they provide some important constraints on stability that must be fulfilled. 

The equatorial thermal and dynamical structure of planetary atmospheres is also 
perturbed by quasi-periodic oscillations. A temperature oscillation with a period of 4-5 Earth 
years, known as the Quasi-Quadrennial Oscillation (QQO), was discovered in the equatorial 
stratosphere of Jupiter 12 , A recent study 13 suggests that there may be a wind oscillation on Jupiter 
that is related to the temperature oscillation, similar to the Quasi-Biennial Oscillation (QBO) on 
Earth. Cassini observations 4,5 recently showed a Semi-Annual Oscillation (SAO) of stratospheric 
temperature in the equatorial region of Saturn. However, we have little information about the 
temporal variation of equatorial stratospheric winds. Such observations would not only enrich 
our knowledge of giant-planet dynamics but also help us explore the full phenomenology of 
Saturn’s SAO. 

Here we use long-term (2004-2009) observations from ISS and CIRS onboard Cassini to 
examine the temporal variations of the equatorial jet in both the troposphere and stratosphere of 
Saturn. In order to obtain accurate measurements of zonal winds in the equatorial region, we 
select only the ISS image pairs with relatively long time separations and high spatial resolutions 
to decrease the uncertainty of wind measurements. In addition, the selected ISS images contain 
long segments of the planetary limb so that they can be navigated well. The selected images on 
two dates (May 10, 2004 and December 30, 2008) and the image processing steps are described 
in the Supplementary Information. Figure 1 shows examples of the processed maps at the strong 
methane filter (MT3) and the corresponding continuum filter (CB3). The MT3 observations are 
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restricted to the hydrocarbon hazes or other photolysis products in the upper atmosphere ' , 
which are ~ 60 mbar in the equatorial region of Saturn 15 . The CB3 observations are generally 
thought to represent ammonia clouds around 500 mbar. The map scale of Fig. 1 is 0.1 “/pixel in 
latitude and longitude, which corresponds to a spatial resolution of - 105 km/pixel at the equator. 
The region very close to the equator (0-5°S) has low contrast due to the scattered sunlight from 
Saturn’s rings. We identified few cloud features in this region, and exclude it from this study. 

After selecting the ISS images, we searched for simultaneous thermal observations in the 
CIRS data. We found CIRS observations on December 29, 2008, nearly simultaneous to the ISS 
observations on December 30, 2008. Unfortunately, we did not find CIRS observations 
corresponding to the ISS observations on May 10, 2004. The closest CIRS observations of the 
equatorial region were taken on April 8, 2005. Details of the selected CIRS observations and the 
procedure of temperature retrieval are also discussed in the Supplementary Information. 

Based on the ISS and CIRS observations, we apply three methods to measure the zonal 
winds in the equatorial region of Saturn. The first method is a manual feature tracking. By 
tracking cloud features that do not substantially alter their shape over time, we can convert their 
movements into the zonal winds. The derived winds are averaged in 1° latitude bins. The second 
method is similar to that of Limaye 16 , in which east-west strips of pixels are extracted from 
image pairs and shifted to find the maximal cross-correlation coefficient. The third method 
applies the thermal wind equation to the CIRS temperature maps, as described in our previous 
study 9 . We use winds measured from the ISS CB3 images to define the boundary level of the 
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thermal wind integration, in order to maintain approximate consistence in time between the 
boundary condition and the CIRS temperature. 

Wind measurements in Saturn’s upper troposphere (i.e., ~50-500 mbar) are shown in Fig. 
2. The uncertainty of wind measurements based on the ISS images comes from three sources: the 
standard deviation of multiple wind measurements, the optical characteristics of ISS cameras, 
and the accuracy of locating cloud/haze features in the ISS images. Uncertainty is discussed in 
the Supplementary Information. Panel A of Fig.2 shows that the two image-based methods give a 
consistent result from the ISS CB3 images. The CB3 zonal wind near 500 mbar increased a few 
m/s from 2004 to 2008, which is within the uncertainty of the measurement. Measurements on 
the ISS MT3 images (panel B), which probe the relatively high hazes around 60-mbar, reveal an 
increase of ~ 20 m/s from 2004 to 2008 around 12°S, significantly larger than the uncertainty. 
We conclude that the wind speed at the level of the high-altitude hazes (~ 60 mbar) increased, 
while the wind speed at the level of ammonia clouds (~ 500 mbar) remained constant (within 
uncertainty) from 2004 to 2008. 

To examine the temporal variation of zonal winds revealed by the ISS images, we further 
analyze the thermal winds in the equatorial region of Saturn. Panel B of Fig.2 shows the profile 
of thermally derived zonal winds at 60 mbar, which is based on Fig. SI in the Supplementary 
Information. The 60-mbar thermal winds also show an increase of ~ 20 m/s around 12°S between 
2005 and 2008. The increase is mainly due to the warming of equatorial region and the 
corresponding changes of the temperature gradient in the troposphere of Saturn (Fig. SI). The 
uncertainty of thermal winds is closely related to the uncertainty of the retri eved temperature and 
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the corresponding meridional gradient, which is difficult to estimate quantitatively. The 
uncertainty of the CIRS retrieved temperature is ~ 1 K in the upper troposphere 17 , smaller than 
the magnitude of temperature variation of ~ 2 K in the equatorial region shown in Fig. SI. The 
retrieval process leaves some uncertainty in the magnitude of warming around the equator, but 
the warming trend from 2005 to 2008 is robust 17 . The warming trend of atmospheric temperature 
intensifies the zonal winds via the thermal wind relationship, which explains the qualitative 
consistency of the high-altitude winds between the ISS measurements and the CIRS observations. 
On the other hand, the temperature of the relatively deep atmosphere around 500 mbar does not 
change significantly (Fig. SI), which is probably due to the long time scale of radiative 
adjustment of the deep atmosphere. 

We also apply the thermal wind equation to the stratospheric temperature retrieved from 
the CIRS nadir observations (Fig. S2 in the Supplementary Information). Fig. S2 shows that the 
thermal winds in the stratosphere increased ~ 60 m/s from 2005 to 2008, which confirms the 
previous results from the CIRS limb observations 18 . The temporal variation of the equatorial jet 
in the stratosphere is related to the varying temperature, which has been observed previously 17 ' 19 . 
In these previous studies 18,19 , the varying temperature and the related thermal winds are further 
linked to the vertical propagation of Saturn’s equatorial oscillation. 

In this study, we report the temporal variation of equatorial jet speed at different altitudes 
in the troposphere and stratosphere of Saturn. We find that the wind speed of the equatorial jet 
around the tropopause (~ 60 mbar) increased ~ 20 m/s from 2004 to 2008. This is the first change 
of the equatorial jet around the tropopause to be observed at the time scale of a few Earth years. 
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At the time scale of one Saturn year (~ 29.4 Earth years), inter-annual variations of cloud 
morphology 20 and heat budget 21 have been discovered on the giant planets. The season of Saturn 
will be the same in 2010-2011 as it was during the Voyager epoch (1980-1981). A detailed 
comparison of the wind field between the Cassini extended mission and the Voyager mission 
will shed light on the inter-annual variation of equatorial jet structure on Saturn. 

As well as enriching our knowledge of the equatorial jet and oscillation, the temporal 
variation of zonal winds has important implications for atmospheric dynamics on Saturn. The 
high-altitude equatorial jet is at 60-mbar, which is near the tropopause 17 . The direct measurement 
of zonal winds near the tropopause will benefit the exploration of thermal winds in stratosphere. 
Fouchet el al. 4 and Guerlet et al. 18 integrated the thermal wind equation in Saturn’s low-latitude 
stratosphere by assuming a null wind as the boundary condition around 20 mbar. Obviously, 
measured winds near the tropopause will be a better boundary condition for deriving the thermal 
winds in the stratosphere, even though the different boundary conditions do not change the 
vertical structure of zonal winds discussed in the previous studies 4,18 . The weaker time variation 
of the equatorial jet at the 60-mbar (~ 20 m/s) than at the 1-mbar (~ 60 m/s) further implies a 
commonality in the vertical propagation of the equatorial oscillation between Saturn and Earth 22 . 
Finally, the difference in behavior in the time variation of zonal winds between the troposphere 
(e.g., 500-mbar) and the stratosphere (e.g., 1-mbar) suggests that the driving mechanisms of 
zonal winds are probably different between the stratosphere, where the radiative effects are 
important and the troposphere where the convection and turbulence are important. 
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Figure Captions 


Figure 1. Southern equatorial maps based on 1SS multi-filter images taken in May 10, 2004. The 
spatial resolution of the cylindrical maps is ~ 105 km/pixel (please refer to Table SI in the 
Supplementary Information). (A) MT3 map at time 1 (tl). (B) MT3 map at time 2 (t2), which is 
taken by ISS in 9.5 hours later than tl . (C) and (D) are the same as (A) and (B) except for the ISS 
images taken at filter CB3. 

Figure 2. Saturn wind measurements on two different dates. (A) Low-altitude winds from the 
ISS/CB3 images. The CB3 winds are located at or deeper than the 500-mbar pressure level, and 
used as the boundary condition of thermal winds. (B) High-altitude winds from the ISS/MT3 
images and the CIRS thermal observations. The CIRS thermal winds at 60 mbar are taken from 
the tropospheric thermal winds in the Supplementary Information. The horizontal solid lines and 
dashed lines are the uncertainties for the cloud-tracking winds and the correlation winds, 
respectively, 
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